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Abstract. We show that the field dependence of the magnetic penetration depth (A), 
for which muon spin rotation (/^SR) is an excellent microscopic probe, provides useful 
information on the degree of anisotropy for the superconducting order parameter. In 
type II superconductors associated with anisotropic order parameters, A is sensitive to 
the quasiparticle excitation induced by the Doppler shift due to a supercurrent around 
magnetic vortices. The presence of such low-energy excitations manifests itself in the 
non-zero slope of A against an external magnetic field. We review recent results on the 
field dependence of A obtained from the application of /xSR to novel superconductors 
that exhibit unconventional characters associated with the anisotropic order parameter. 



PACS numbers: 74.25.Ha, 74.70.-b, 76. 75. -hi 
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1. Introduction 

One of the greatest impacts from the discovery of high-Tc cuprate superconductors 
is renewed interest in the exotic mechanisms of superconductivity, and a boosted 
search for their model systems, which has led to a soaring number of novel materials 
identified as superconductors [T]. They include varieties of transition-metal oxides, 
borides, borocarbides, and other intermetallic compounds with rare-earth elements. 
The list can be readily extended by including those based on organic materials, and 
is still growing in length. Compared with classical simple metals or binary compound 
superconductors, they have a common distinctive feature that the pairing correlation is 
potentially highly anisotropic due to a strong electronic correlation (Coulomb repulsion) 
and/or a considerably two-dimensional nature of the Fermi surface, which is also shared 
by cuprate superconductors. 

It is well established that any metallic system can fall into a superconducting state 
when there is an effective interaction that is attractive between the conduction electrons 
(pair correlation). In this situation, electrons tend to form a bound (pairing) state 
between those with opposite momenta (k and — k), so that the phase of the pairing 
wave-function may become absolute zeroj2]. This instability towards the formation of 
electronic bound states without a barrier is an intrinsic character of the Fermi surface, 
irrespective of the microscopic mechanism of the attractive interaction. When the 
pair correlation is mediated by the electron-phonon interaction, as in ordinary cubic 
metals, the pair correlation has the least dependence on k (relative orbital angular 
momentum L = 0), and thereby the structure of the superconducting order parameter 
A(k) is isotropic over the entire Fermi surface. Because of the mandatory requirement of 
Fermi statistics that the electronic wave-function must be anti-symmetric, the remaining 
freedom of the spin state in the paring function is set to be a singlet {S = 0). This pair 
correlation, having s-wave and spin singlet symmetry, is conventionally called the "BCS 
mechanism"; in its narrower sense, the electron-phonon interaction is presumed to be 
the primary origin of the attractive interaction [2]. 

Historically, the first sign of non-BCS type pairing was found in the superfluidity 
of liquid ^He^ Ej, where the neutral ^He atoms (which have a nuclear spin of 1/2 
and thereby obey Fermi statistics) play the role of electrons in superconductivity. It 
is known that liquid '^He can be regarded as being a Fermi liquid (i.e., having a well- 
defined Fermi surface) below ~ 10^ mK, and that it becomes a superfiuid below a few 
mK, where it can flow narrow channels without friction (superfluidity). However, there 
is a major difference between the nature of pairing between ^He atoms and that of the 
BCS type. The ^He atoms have a hard core of a repulsive interaction with a relatively 
large radius, which makes it difficult to pair in a state with zero angular momentum. 
Thus, many theories have predicted that ^He atoms may pair in a p-wave (L = 1) or d- 
wave {L = 2) state, where they can keep themselves apart while the pairing interaction 
is at workjniEI. Later, experiments confirmed that they were indeed in a p-wave state. 
This also meant that the spin part of the pairing wave-function (order parameter) must 
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be in a triplet state {S = 1), leading to a variety of possibilities for the total state of the 
pairing wave-function (oc A(k)) to break symmetry. It is now established that there are 
at least three different phases in the superfluidity of liquid ^He that all have different 
order parameters. 

It is readily predicted that the situation similar to that in liquid ^He can be 
realized when a short-range repulsive electronic correlation is not negligible in metallic 
superconductivity. As has been established during the past decade, high-Tc cuprates are 
among the first such examples in which electrons pair in a state other than an s-wave 
due to a strong on-site repulsive correlation. The Cooper pairs in cuprates prefer a 
(i-wave because of the tetragonal structure of the two-dimensional Cu02 lattice and the 
associated symmetry of the Fermi surface. More specifically, the pairing is a (i^2_y2-wave 
and the order parameter is described by a gap function. 



where Aq is the maximum value of the anisotropic gap; the energy gap vanishes along 
the lines kx = ^ky, which are called line nodes. It also happens that the electronic 
correlation in cuprates is antiferromagnetic, as is naturally expected for doped Mott 
insulators, which makes it favorable to form spin singlet pairs. The latter points to a 
magnetic origin as a pairing mechanism, irrespective of the true nature of the ground 
state which is still under debate. Thus, the structure of the superconducting order 
parameter reflects important characteristics of the electronic correlation. 

In this paper, we review our recent studies on the structure of the superconducting 
order parameter in novel type II superconductors by muon spin rotation/relaxation 
(/iSR). It is well known that a magnetic field can penetrate type II superconductors as a 
bundle of quantum flux hues (magnetic vortices) jH] , where the spatial field distribution, 
B[r), becomes inhomogeneous due to gradual change in the supercurrent flow around the 
vortices. The degree of inhomogeneity, which is primarily determined by the magnetic 
penetration depth (A), a magnetic cutoff parameter (^v, which is proportional to the 
Ginsburg-Landau (GL) coherence length ^ql) and the spacing of vortices (a) can be 
measured directly by /xSR as a spin-spin relaxation (I/T2). By applying a refined 
analysis technique, one can reconstruct B{r) more accurately so that both A and .^v 
may be deduced separately jH]. Among various experimental techniques applied to a 
similar end, /iSR technique is unique in many respects. Namely, it is a microscopic 
technique that can be applied to virtually any superconductors having a reasonable 
magnetic penetration depth (A < 5000 A). Because of the purely magnetic nature of 
a spin 1/2 probe, the interpretation of the /xSR spectra is free from any complication 
due to additional interactions from higher multipoles often found in nuclear magnetic 
resonance (NMR). The /iSR technique is sensitive to the bulk property, and is thus free 
from effects specific to the surface, while they are often problems for scanning tunneling 
spectroscopy (STS) or angular- resolved photo-emission spectroscopy (ARPES). 

In the following, we demonstrate that the temperature/field dependence of A is 
strongly affected by the anisotropy of the order parameter. In particular, A is enhanced 
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by an external field due to the Doppler shift of quasiparticles in the gap nodesjTU], which 
leads to almost a linear increase of A with increasing field. This feature is regarded as 
an unambiguous sign of the presence of nodes in the energy gap. We show several 
examples of field-dependent A, some of which are indeed identified to have anisotropic 
order parameter by other experimental techniques. A more comprehensive review of a 
similar study can be found elsewhere jH]. 

2. Internal Magnetic Field Distribution in the Mixed State 
2.1. Microscopic Model 

In penetration-depth measurements, it is common to assume a geometrical condition 
that muons are implanted into a specimen with the initial spin polarization 
perpendicular to the external field, H = {0,0, H). Then, since muons stop randomly 
along the length scale of the flux line lattice (FLL), the time evolution of complex 
muon polarization, P(t), provides a random sampling of the internal field distribution, 
B(r) = (0,0,5(r)): 

P{t) = P,{t) + iPy{t) = e^-aie) / n{B) expit^^Bt + (l))dB, (2) 

J —oo 

n{B) = {5{B - B{v)))^, (3) 

where Px,y{t) is proportional to the time-dependent ii^-e'^ decay asymmetry, A^^yif), 
deduced from a corresponding sets of positron counters, n{B) is the spectral density 
for the internal field defined as a spatial average (( )r) of the delta function, 7^ is 
the muon gyromagnetic ratio (=27rxl35.53 MHz/T), and (p is the initial phase of 
muon precession im IT^ . The additional relaxation (ab) is mainly due to random local 
fields from nuclear magnetic moments (un ~0.1 yus~^) and the distortion of flux line 
lattice due to a random pinning of vortices (cTp), which can be approximated by a 
Gaussian relaxation^Jj. These equations indicate that the real amplitude of the Fourier- 
transformed muon precession signal corresponds to n{B) with an appropriate correction 
of (Tb. While (T„ can be estimated from the spectrum in the normal state, dp often needs 
to be considered as a variable parameter, depending on the temperature and/or external 
magnetic field. 

In modeling the internal field distribution of vortex state, the simplest approach is to 
assume that the field distribution is a linear superposition of that for an isolated vortex, 
as presumed in the London theory |T3]. This is a reasonable assumption as long as the 
inter- vortex distance is much longer than the GL coherence length {H <^ -f^c2)- Then, 
what we need is to know the field distribution around a single vortex and the structure 
of the vortex lattice. The latter must be given from other sources of information, such as 
small angle neutron scattering (SANS) or scanning tunneling microscopy/spectroscopy 
(STM/STS). In the London model, -B(r) is approximated as the sum of the magnetic 
induction from isolated vortices, to yield 

^(r) = i?oEY^^^(^'^v), (4) 
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Figure 1. A magnetic field distribution n{B) calculated by the modified London 
model with a triangular flux line lattice. (The Lorentzian cutoff was assumed.) The 
marking 'c' refers to the contribution from vortex cores, 's' to saddles, and 'v' to valleys, 
respectively. The Gaussian distribution is only meant for a guide to eye. 



where K are the vortex reciprocal lattice vectors, Bq (~ H) is the average internal field, 
A is the London penetration depth, and F{K, C,v) is a nonlocal correction term with 
being the cutoff parameter for the magnetic field distribution; care must be taken not 
to interpret naively as ^gl which is for the spatial variation of the superconducting 
order parameter. While the Lorentzian cutoff, F(K,^^) = exp(— a/2K^v), is predicted 
to be a better approximation for the GL theory at lower fields^!], the Gaussian cutoff, 
^v) = exp(— generally provides satisfactory agreement with the data[T^. 
Note, however, that the Gaussian cutoff is derived from the GL equations near Hc2, 
and thus would not be appropriate at lower fields. A comparison of the analysis on 
identical data indicates that a Gaussian cutoff yields a signficantly larger value for A 
and a stronger field dependence (about a factor 2) than those obtained by a Lorentzian 
cutoff p. 

Besides the London model, there are a couple of models for B{r) based on the 
Ginsburg-Landau theory. Although the GL equations can be solved to yield an 
approximate analytical solution for the mixed state near Hd or Hc2, it must be solved 
numerically for intermediate fields. Fortunately, it is known that the field distribution 
obtained from exact numerical solutions of the GL equations is in excellent agreement 
with that from the modified London model at low fields and arbitrary k (= A/^gl, the 
GL parameter) [THj . 

One of the most important characteristics of n{B) obtained from these models is 
the site-selective feature of the line shape. As shown in Fig. the sharp peak due to the 
van Hove singularity found in the lower field mainly represents the contribution from 
the saddle points of B{r), the lower field end from the central valleys among vortices. 
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and the high field end from the region near the vortex cores. This quite asymmetric field 
profile with such a geometrical correspondence allows us to determine A and reliably 
by comparing the time evolution of the muon spin polarization with that calculated by 
n{B). 

It is often the case that the superconducting properties in non-cubic compounds 
are strongly anisotropic, leading to a large difference between the magnetic penetration 
depths for in-plane and perpendicular directions. More specifically, in uniaxial 
superconductors with M^h and Mc being the carrier mass for in-plane and perpendicular 
directions, we have 

B{r, e) = BoY: KK)e-^^-^F(ir, ^v), (5) 

K 

^ ^ (1 + K2mabA2 + ir>..A2) (1 + K^m^A^) - K^K^ml^X^ ' ^ ' 

where 6 is the polar angle of the c-axis, = Mg}^/M, rric = M^/M (with M = 
(M^MJVS), and 

ITT'xx = ITT'ab COS^ 6 + nic siu^ 9, (7) 

^zz = ^ab sin^ 9 + nic cos^ 9, (8) 

"fTT-xz = {f^ah — ^c) Sm9cOs9. (9) 

Thus, the line shape depends on Aab = Ay^m^ and Ac = \^/rn^ in a complex manner [JTj. 



2.2. Gaussian field distribution 

When the quality of /xSR data is good enough to be analyzed by the above model, we 
can obtain A and simultaneously by directly comparing the //SR time spectrum with 
that calculated from -B(r). Unfortunately, our experience shows that this is not always 
the case when the sample is not a single crystal, or A happens to be very large, etc. so 
that the characteristic features of n{B) and the associated time spectra important for 
such analysis are smeared out. In such a situation, the Gaussian field distribution has 
been used as a convenient analytical model, where the depolarization rate is presumed 
to be given by the second moment of the field distribution (A 3> ^v), 

{^B') = mv) - H)\, (10) 

which is reflected as T2 relaxation in the /iSR line shape. The Gaussian distribution of 
local fields naturally leads to a Gaussian depolarization function, 

P{t) ~ exp(-(T^t2) exp(-(72tV2) exp{i-f^Ht + 0), (11) 



a =^,^{AB^). (12) 

For the ideal case of a triangular FLL with an isotropic effective carrier mass and a cutoff 
K ^ 1.4/^v provided by the numerical solution of the GL theory, A can be deduced from 
cr using the following relation [THl IT^ IT^ : 

a{h) [fxs-^] = 4.83 X 10^(1 - h)\-^ [nm], (13) 
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Figure 2. Magnetic penetration depth (Ag) deduced from an analysis using the 
Gaussian damping (fiUed triangles and circles), where the analyzed spectra were 
those simulated by the modified London model with A shown by the open triangles 
(corresponding to 77 = 0) and circles (77 — 2) at each field. A linear fit yields ry = 0.11 
and 1.53, respectively. The upper-critical field was Hc2 — 10 T, from which the cutoff 
parameter (^v) was determined to be 40 A ( = 0-7^gl, where ^ql = \/ ^q/2'kHc2, with 
$0 being the quantum flux). 



where h = H/Hc2- While the above form is vahd for h < 0.25 or h > 0.7, a more useful 
approximation valid for an arbitrary field is ^2] 

a{h) [i2s-^] = 4.83 X 10^(1 - h)[l + 3.9(1 - hy]^/^X-^ [nm]. (14) 

The field dependence of those equations represent a reduction of the Gaussian width 
due to a stronger overlap of vortices at higher fields, while A is a constant; therefore, 
the deviation of a{h) from those equations can be attributed to the change of A with 
the field. 

However, the microscopic situation of the FLL state is considerably different from 
the above "ideal" one in practical cases where Gaussian damping is actually observed, 
because there must be an additional effect of randomness to round up the sharp feature 
of n{B). This also makes it difficult to distinguish ab from a in Eq. ()11|). giving rise to 
a problem in comparing the values of A, for example, between those from an analysis 
using the modified London model and those obtained from the Gaussian approximation 
(where the influence of cJb is indistinguishable). 

In order to examine the model dependence of the analysis, we made a simulation 
to compare the results of an analysis, where we generated /iSR time spectra using a 
modified-London model, and then analyzed them by simple Gaussian damping (Eq. |TT|) ) 



Field-induced quasiparticle excitations in novel type II superconductors 



8 



to deduce a. According to our result, the Gaussian distribution originates from the 
distribution of A (which may vary at different sample domains), which must be present 
in the polycrystalline powder specimen of anisotropic superconductors (e.g., Ac ^ Aab)- 
As is also clear in the field profile shown in Fig.Q the simulated time spectra with typical 
values for A (2000-3000 A) cannot be fitted by Eq. (fTT|) due to the strongly exponential- 
like damping; this is obviously due to the contribution of high-frequency tails in the 
spectral density, n{B). (Thus, the use of the second moment as an approximation in 
the ideal situation would be valid only when the relaxation rate is small enough to 
eliminate the asymmetric feature of n{B).) The situation was much improved when the 
Gaussian distribution of A was introduced with a variance {(Tx), 

G(A)(xexp[-(A-A)Va^], (15) 

where A is the mean value. For the parameter values shown in Fig. |21 the time spectra 
become Gaussian-like when ax ~600-800 A, yielding reasonable reduced chi-square 
values by Eq. (fTT|) . We also assumed a gradual decrease of ax, 

axih) = axil-h''), (16) 

considering that the elastic moduli Ca of FLL, which control the FLL distortion and 
the associated modulation of A, exhibit a quadratic dependence on the applied field. 
(Note, however, that the factor 1 — h? yields only a small change of ax{h) for h < 0.5.) 
Some examples are shown in Fig. |21for ax = 800 A, where the penetration depth (Aq), 
obtained by Eq. (fT^ from a{h), is plotted together with the original A. A reasonable 
agreement between Ag and A is seen, except for the case when A = 2000 A, where 
Ag takes systematically lower values at all fields. This can be readily understood by 
considering the fact that a is enhanced by an amount dp as a consequence of Eq. (fTHjl . 
We also examined the field dependence of A, which would be most crucial in the following 
arguments, 

Xih) = XiO)[l + r]-h], (17) 

where ?7 is a dimensionless parameter used to express the magnitude of the field 
dependence. As is evident in Fig. |21 the slope dXc/dh is slightly weaker than the 
original assumption; when we take t] = 2 for the simulation, we obtain t] = 1.53 as the 
corresponding slope for Ag. However, when there is no field dependence of A (r/ = 0), Ag 
exhibits the least dependence on the field (?7 = 0.11). Thus, we can conclude that the 
field dependence of the penetration depth (as a mean value) obtained from the Gaussian 
field approximation provides a sound basis for the characterization of superconductors. 

3. The Doppler shift and associated non-linear effect 

In the FLL state, the quasiparticle momentum vp is shifted by the flow of supercurrent 
Vg around the vortices due to a semi-classical Doppler shift, leading to a shift of the 
quasiparticle energy spectrum to an amount e = mv-p ■ Vg. Since the density of state 
(DOS), N{E), is non-zero, except at the Fermi level {E = 0) and is higher at larger 
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Figure 3. Fermi surface shifted by the quasiclassical Doppler effect due to the 
supercurrent around vortices. While such a shift has no effect on the quasiparticle 
excitation for the isotropic gap (£(vf • Vg) < As), it induces additional excitation by 
breaking pairs near to the nodal region (pointed by arrows). 



energy {0 < E < Aq) for the anisotropic order parameter, quasiparticles can be excited 
by the Doppler shift outside of the vortex cores with a population proportional to 
N[E + e{vp ■ Vs)), leading to an enhancement of A [lOj. In other words, the Cooper 
pairs with a gap energy of less than e can be broken by the Doppler shift (see Fig. E}. 
Historically, a similar effect was considered first for type I superconductors, where the 
non-linear response of the shielding current in the Meissner effect due to the "backflow" 
of quasiparticles was discussed pUj. 

The magnitude of rj represents the degree for the increase of DOS for quasiparticles, 
which must be roughly proportional to the phase volume of the Fermi surface where the 
Doppler shift exceeds the gap energy (£:(vf ■ Vg) > A(k)). It also follows that the effect 
depends on the direction of Vg (and hence that of the external field H relative to the 
order parameter) in a sigle crystalline specimen. According to Volovik, the quasiparticle 
density of state for the anisotropic order parameter is 

iVdeioc(O) - NpKecL^ = Npgih), (18) 
K^f |A(k)Mk, (19) 

where Np is the DOS for the normal state and i^' is a constant on the order of unity ^Oj- 
It is important to note that K is proportional to the phase volume of the low excitation 
energy in A(k), thereby carrying information on the degree of anisotropy for A(k); the 
factor /i^/^ comes from the inter-vortex distance (oc h~^/'^) multiplied by the number of 
vortices (oc h). The superfluid density at a given field is then 

ns(/i) ~n,(0)[l-(7(/i)], (20) 




Figure 4. Field dependence of A normalized by the value at /i = for the case 
g{h) = y/h (A'Cgl = !)■ It is well represented by a linear dependence, as in Eq. H22II 
with c ~ 1.5. 



which is directly reflected in the magnetic penetration depth, 
1 47re2 



risih). (21) 



Therefore, as a mean approximation, we have 

m = -^== - A(0)[1 + cKeo^hl (22) 

where c ^ 1.5 ioi < h < 0.5, as show in in Fig. IH Thus, the comparison between 
Eqs. dni) and ^ yields 

V ^ cK^GL, (23) 
indicating that the slope t] reflects the phase volume of the Fermi surface where 

|A(k)|<5. 

Since the Doppler shift is far smaller than the gap energy in the usual situation 
for the isotropic gap, no such enhancement is expected for the conventional s-wave 
pairing [t] <^1). A recent theoretical calculation based on the Bogoliubov-de-Gennes 
(BdG) equations indicates, however, i] is not exactly zero for s-wave pairing, although 
it is much smaller than that for (i^2_.^2 -pairing pi] . It must also be noted that the 
effect of temperature must be considered to evaluate the degree of anisotropy from the 
measurement of rj. In general, the phase volume of the Fermi surface where |A(k)| < ksT 
also contributes to quasiparticle excitation. Thus, strictly speaking, the observation of 
a finite rj means the presence of a small gap region with an upper bound of ~ ksT (see 
below, the case of MgB2 for example). 
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4. Non-local corrections 

Because a superconducting pair correlation occurs over a finite length scale, (i-e, the 
BCS coherence length), the electromagnetic response of superconductors is subject to 
various non-local effects. The primary example is the cutoff term, F{K, C,v), incorporated 
in the modified London model (Eq. (^). Moreover, there are a couple of other 
corrections that must be considered for anisotropic superconductors. 

In the superconducting state with gap nodes in the order parameter, the 
quasiparticles are mostly confined to the vicinity of the nodes at low temperatures. 
This generally leads to the suppression of quasiparticle excitation due to the non-local 
electrodynamics caused by the divergent coherence length, ^ oc A(k)~^, yielding a 
weaker temperature/field dependence of A at higher fields f22[ I23j. Such an effect has 
been studied experimentally in detail for the case of YBa^CuciOK or^j^ IH]. 

Another important correction comes from the anisotropy of the Fermi surface. 
Novel superconductors including cuprates have a common feature that the Fermi 
surface tends to exhibit a strong anisotropy due to, e.g., a two-dimensional and/or 
multi-band character, which influences the flow of supercurrent over a length scale, 
^VF/Aoj2Sl I2ni I2H]- More speciflcally, the length scale is also controlled by the 
mean-free path (/) for electrons. It turns out that this non-local correction partially 
accounts for the change in the vortex lattice structure from triangular to squared lattice 
in various systems, including i^NigBaC (/2=Y, Lu) 123 EDI ED E3 El , VgSilSaEni, 
and La2-xSr2:Cu04 (x = 0.17)[37 . This also leads to a change in -B(r) due to the 
modifled flow of supercurrent from circular to squared shape (see below, e.g. YNi2B2C). 
Note, however, that the difference in the free energy between a triangular and a 
square vortex lattice is fairly small, making the lowest-energy conflguration strongly 
dependent on other physical quantities, such as the temperature, magnetic fleld, and 
crystal orientation, which is also in a strong correlation with the superconducting order 
parameter. 

5. Overview of /iSR results 

In this section, we try to establish the correspondence between the presence of the 
anisotropic order parameter and a non-zero slope (rj) in the magnetic-fleld dependence 
of X{h). As shown below, rj provides a good measure for the degree of anisotropy in the 
superconducting order parameter. However, one has to keep in mind that additional 
information is generally needed to resolve the precise symmetry of the pairing; one would 
be easily led to a false conclusion in choosing, e.g., between d- and s + g-wave pairing 
based solely upon the /xSR result. Another origin of apparent anisotropy would be 
multi-gap order parameters with one of those having a small gap (< k-sT), as suggested 
in the case of MgB2 (see below). 

While a vast body of superconductors have been investigated by /xSR, there are not 
many of them in which the fleld dependence of X{h) has been measured in detail. This 
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Figure 5. Field dependence of the magnetic penetration depth (A) in CeRu2 at 2 K 
obtained by fitting data with the modified London model, where the dashed line is a 
guide for the eye (after Ref.03|)- 



is partly due to the historical reason that the /iSR apparatus with a high magnetic field 
has become available for routine service only since the late '90s. Here, we summarize 
our recent work on CeRu2, Y(Pt,Ni)2B2C, Cd2Re207, and MgB2, in which detailed /xSR 
measurements have been performed. The results of NbSe2 and YBa2Cu306.95 are also 
mentioned for a comparison. A wider variety of compounds in view of other experimental 
techniques are covered elsewhere ^35] . 

5.1. CeRu2 

The cubic Laves phase compound, CeRu2 (Tc = 6.1-6.5 K at zero field, Hc2{T = 0) ^ 6- 
7 T), has a long history of experimental and theoretical studies since its discovery in 
1950s. One of the current issues is its magnetic response at higher fields [h > 0.5), 
where an anomalous enhancement of quasiparticle excitation has been reported. This 
is suggested by the observation of de Haas-van Alphen (dHvA) oscillation over a field 
region where the cyclotron radius is much larger than the inter- vortex distance [SHI- The 
presence of excess quasipaticles has been further confirmed by a strong enhancement of 
A measuremed by /iSR at higher fieldsllO]. Moreover, while most of experimental studies 
concluded that the pairing symmetry is a spin-singlet s-wave, detailed studies on the 
spin-lattice relaxation in nuclear quadrupole resonance (NQR) suggested the presence 
of anisotropy in the order parameter |41j. This was apparently in line with the observed 
non-linear field dependence (oc /i^/^) of the electronic specific heat coefficient '7(/t)|42j: 
As indicated in Eq. ([T^ . the quasiparticle excitation has a contribution proportional 
to /i^/^ for the anisotropic order parameter, while an /i-linear dependence is expected 
for the conventional case, because 7(/i) must be proportional to the volume of normal 
cores, and thereby to the number of vortices. 
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However, our /iSR studies in a high quahty single crystal has shown that the field 
dependence of •y^h) can be attributed to that of the vortex core radius Pv (oc .^v), at 
least over the region h < 0.5; the fiSR spectra were analyzed (in the frequency domain) 
by the modified London model to extract and A independently (We later re- 
analyzed a part of the data in time domain and found that the result was unchanged.) 
When the quasiparticles are confined within the vortex cores, the electronic specific- 
heat coefficient must be proportional to the cross section of the cores multiplied by the 
number of vortices. 



where our result for Pv{h) is well reproduced by /3 ~ 0.53. Thus, the observed nearly 
/i^/^ dependence of 7(/?-)|121 exhibits perfect agreement with our result. This is also 
quahtatively in line with a recent theoretical calculation for s-wave superconductors 
based on the quasi-classical Eilenberger equations, where /? ~ 0.67 is predicted due 
to various non-trivial effects, including those from the inter-vortex interaction and the 
vortex core excitation at finite temperature This strongly suggests that quasiparticle 
excitation is confined within the normal cores of the vortices. More importantly, we 
found that X{h) exhibits the least dependence on h (namely, 77 ~ 0) over the relevant field 
range, as clearly shown in Fig. I5(43j . This is perfectly in line with the above conclusion 
obtained for the vortex cores, as is also the case with s-wave pairing suggested by the 
vast majority of other experimental results. 

Meanwhile, it must be stressed that the situation changes drastically upon 
increasing the field above h ~ 0.5. Figure El indicates that A exhibits a divergent 
increase for h > 0.6 with increasing field, thereby suggesting a divergent increase of 
quasiparticle excitation. This is strongly supported by the observation of the dHvA 
effect, as mentioned above. To our knowledge, there is no simple explanation of such 
an anomaly j45j. One possible model may be that proposed by Fulde-Ferrel-Larkin- 
Ovchinnikov, where a new superconducting phase with a spatially inhomogeneous 
order parameter is predicted to occur in rare-earth compounds having a large spin 

paramagnetism |3ni SZ| • 

Finally, we note that the absence of a clear coherence peak in NQR may be due 
to the weak random magnetism observed below ~40 K by zero field /iSR[48 . They 
observed an increase in the muon spin relaxation rate on the order of 0.02/is~^ in 
accordance with the increase of the ac-susceptibilityjini- The sample quality suggests 
that the weak magnetism is of intrinsic origin, which would act as a scattering source 
for pair breaking. 



The borocarbide superconductor, YNi2B2C (Tc = 15.4 K at zero field, Hc2(T = 0) ^ 7-8 
T), has attracted much attention due to the strong /i^/^ dependence of the electronic 
specific-heat coefficient in a high-purity specimenjSH], with which they suggested a 



7(/i) oc Tvplih) ■ h = irplh'^, 

p^h) = poh^^-'y\ 



(24) 
(25) 



5.2. Y(m,Pt)2B2C 
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change of the vortex core radius based on earher experimental suggestions that an 
s-wave pairing was reahzed in this compound; the system showed httle sensitivity to 
non-magnetic impurities, as is typically found in BCS s-wave superconductors. They 
also found that such a /i^/^ dependence was replaced by a /i-linear dependence upon 
substitution of Ni by Pt (~ 20%). By now, there is mounting evidence that the 
order parameter in pure YNi2B2C is considerably anisotropic 1^ 13^ 15^ 1331 13^ I37j . 
although the pairing symmetry is basically s-wave-like. The key to understand the 
varying results in borocarbides is that the anisotropy of the order parameter is indeed 
sensitive to impurities, which is in good contrast to the robustness of superconductivity, 
itself; s + (^-pairing, for example, changes into an effectively isotropic s-wave pairing, 
where the anisotropic part (associated with gf-component) is washed out by impurity 
scattering. Thus, a part of the divided results obtained by NMR-l/Ti measurements 
may be sorted out in terms of the sample purity 38]. Our /zSR study was quite successful 
to clarify the effect of non-magnetic impurities on the anisotropic order parameter in 
Y(Nii_,Pt,)2B2C. 

As mentioned before, there is a significant contribution of the non-local effect in 
borocarbides due to the anisotropy of the Fermi surface. This effect must be considered 
for modeling of the magnetic field distribution, -B(r), and the associated spectral density, 
n{B), to obtain the correct values for A and ^v- To this end, the London model is further 
modified to yield 

where the terms proportional to K^Ky represent the non-local effect, with q being the 
parameters coming from the band structure [21] ■ Moreover, the non-local effect leads 
to the formation of a squared flux line lattice, which also modifles B{r). Our result 
indicates that these features have a strong influence on n{B) probed by /iSR[^3j. For 
example, no reasonable flt can be obtained when one assumes a square FLL without 
non-local correction terms. 

Our nSR experiment has revealed that the vortex core radius, p^,m a. pure specimen 
[x = 0) exhibits a much steeper decrease with increasing fleld than that estimated 
from the electronic speciflc-heat coefficient, 'y{h); p^Qi) shrinks sharply for h < 0.1, 
then changes only very weakly with the field. This is in marked contrast with the 
case of CeRu2, where the field dependence of 7(/i) is in good accord with Pv{h) (i.e., 
^{h) (X hiiplih)). This suggests the presence of excess quasiparticle excitation outside of 
the vortex cores, which contributes to 7. As shown in Fig.lHl this is indeed supported by 
the observation that A exhibits a strong field dependence in a pure specimen. The slope 
is deduced from the linear fitting (Eq. p7j) ) to yield rj = 0.95(1). On the other hand 
in the Pt-doped specimen [x = 0.2), A is mostly independent of the field {rj ~ 0)|58j. 
This is again consistent with the presumed hybrid nature of the order parameter, where 
the s-wave component is relatively enhanced by impurity scattering. We also note that 
A in the Pt-doped specimen is about 1.27-times longer than that expected solely by 
the conventional impurity effect jHE]- This strongly suggests that there is an excess 
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Figure 6. Field dependence of the magnetic penetration depth (A) in 
Y(Nii_a;Pt2:)2B2C, obtained by fitting data with the modified London model and a 
non-local correction, where the data were obtained at 3 K for a; = 0.0 and 2.5 K for 
X = 0.2 (Tc = 12. f K), respectively (after Ref-ISIEHl)- 

quasiparticle density of states generated by the interaction between the impurities and 
the anisotropic component of the order parameter, as is found in superconductors with 
gap nodes. 

5.3. MgB2 

The revelation of superconductivity in a binary intermetallic compound, MgB2, has 
attracted much interest because it exhibits an almost two-times higher transition 
temperature (Tc ~ 39 K) than those of all intermetallic superconductors known 
to datej^. The most interesting issue concerning this compound is whether or 
not it belongs to the class of the conventional BCS type (namely, phonon-mediated 
spin-singlet s-wave pairing) superconductors. So far, most experimental results 
favor phonon-mediated superconductivity [Hill ED 1021 EHl IMl lEEl EHl EZ|- On the 
other hand, calculations of the band structure and the phonon spectrum predict 
a double energy gap jHEl 1001, with a larger gap attributed to two-dimensional 
Px-y orbitals, and a smaller gap to three-dimensional Pz bonding and antibonding 
orbitals. The experimental results of specific heat measurements [ZOIIZII, point-contact 
spectroscopy j72I, photoemission spectroscopy jT^I, scanning tunneling spectroscopy [71] 
and penetration depth measurements jTHI have supported this scenario. 

The double energy gap would have a direct relevance on the temperature 
dependence of A, because there must be excess quasiparticles excited over the smaller 
energy gap (As) at higher temperatures, while the bulk superconductivity is maintained 
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Figure 7. Field dependence of the magnetic penetration depth (A) in MgB2 obtained 
by fitting data with the Gaussian field distribution model, where the data were obtained 
at 10 K. The values for h < 0.1 are probably underestimated due to addtional relaxation 
induced by flux pinning (after r!,ef.|77jV 

by the larger gap (Al). At this stage, there are two such /iSR measurements reporting 
the result of an analysis based on the two-gap model, where Ag is reported to be 2.6(2) 
meVpB] and 1.2(3) meVjZZj, respectively. On the other hand, the field dependence 
of A is sensitive only to those excited by the Doppler shift, and therefore the slope rj 
would be zero as long as both energy gaps are isotropic. The only exception would 
be that the temperature at which X{h) is measured is comparable to A^/k-Q, so that 
the smaller gap effectively becomes equivalent with nodes in the order parameter. As 
shown in Fig. [3 our result of X{h) exhibits a clear dependence on h with r] = 1.27(29), 
where the measurements were performed at T ~ 10 K^^. Considering that the 
measured temperature nearly corresponds to As/k^ = 14(4) K, the observed field 
dependence would be qualitatively consistent with the two-gap model with isotropic 
order parameters. Meanwhile, if the smaller gap is as large as 2.6 meV|76j. it would 
mean that either Al or Ag is anisotropic. 

Unfortunately, so far it is difficult to obtain a single crystal of MgB2 with the 
dimensions necessary for the conventional /xSR experiment; thus, all of the /xSR 
measurements have been performed on powder specimens. The obtained time spectra 
were fitted by the Gaussian field distribution (Eq. (jllj) ). where the additional relaxation 
due to the fiux pinning (oc exp(— cXpt^)) was not separable. As a matter of fact, we 
observed an increase of the relaxation rate with increasing field over the region h < 0.1 
in our MgB2 specimen, which might be related to the fiux pinning [77]. However, as we 
showed in Section 2.2, an analysis based on the Gaussian field distribution model has a 
relatively weak uncertainty in terms of the relative change of A against the field. Thus, 
we think that the above result provides a sound basis for a qualitative evaluation of the 
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gap anisotropy. 

Another source of complication would be that there is a strong anisotropy of 
magnetic response over the crystal direction, which is mixed up in the polycrystalline 
specimen; it has been reported that the upper critical field for H parallel to the a6-plane 
(11^2) is about three-times higher than that for the perpendicular direction (ij"j^n)|7^. 
This introduces an uncertainty in the definition of the normalized field h (= H/Hc2), 
which is directly reflected in the evaluation of rj. Thus, further measurements on a single 
crystalline specimen would be necessary for the reliable evaluation of rj. 

Recently, we made /iSR measurements on a new superconductor, Ca(Sio.5Alo.5)2 
(Tc = 7.7 K), which has a crystal structure quite similar to that of MgB2[79j. Provided 
that the structure of order parameter in MgB2 is strongly related with that of the Fermi 
surface, a similar situation might be expected in this compound. Our preliminary result 
on a polycrystalline specimen with the Gaussian analysis indicates that A exhibits a 
field dependence with rj ~ 1.85, thereby supporting the above conjecture, at least in 
terms of quasiparticle excitation [HHj. 

5.4. Cd2Re207 

A class of metal oxides isostructural to mineral pyrochlore has been attracting 
considerable attention because they exhibit a wide variety of interesting physical 
properties jHJ. The pyrochlore has a general formula of A2B2O7, consisting of BOg 
octahedra and eightfold coordinated A cations, where A and B are transition metals 
and/or rare-earth elements. In particular, the B sublattice can be viewed as a 
three-dimensional network of corner-sharing tetrahedra, providing a testing ground 
for studying the role of geometrical frustration in systems that have local spins at 
B sites with an antiferromagnetic (AFM) correlation. Although metallic pyrochlores 
comprise a minority subgroup of the pyrochlore family, they consist of distinct members, 
such as Tl2Mn207, which exhibits a colossal magnetoresistance. \S3\ IHl] In view of these 
backgrounds, the recently revealed superconductivity in 5d transition metal pyrochlores 
and related oxides, Cd2Re207|H3Eni and KOs206[S7j, is intriguing, because they evoke 
anticipation for exotic superconductivity. 

It is reported that Cd2Re207 falls into the bulk superconducting state below 
Tc ^ 1 ~ 2 K[S3]. The dc- magnetization curve indicates that the superconductivity 
is of type II with the upper critical field close to 0.29 T at K. So far, ^^^Re NMR-l/Ti 
exhibits a clear coherence peak typically found for the conventional s-wave pairing jHH], 
although there are not many reports concerned with the pairing symmetry. Thus, the 
existing evidence strongly suggests that the order parameter is unexpectedly isotropic. 
This is further supported by the absence of a field dependence for X{h). Figure |H1 shows 
X{h) versus h, where one can clearly see that ~ over the observed field range [HH]. 
Here, we note that the upper-critical field can be determined by a /iSR measurement 
when it is well within the reach of the apparatus (< 7 T). As is evident in Eqs. (jlHj] 
and ()14|) . the spin relaxation due to FLL is quenched at /i = 1 (i.e., H = ifc2)- Thus, 
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Figure 8. Field dependence of tlie magnetic penetration depth (A) in Cd2Re207 
obtained by fitting data with the Gaussian field distribution model, where the data 
were obtained at 0.2 K (after Ref.|^). 

from the field dependence of a, we obtained Hc2 = 0.37(5) T for our specimen. The 
normahzed field in Fig. |H1 is defined by this value for Hc2- The fact that A exhibits 
the least dependence on h also means that the field dependence of a is well reproduced 
by Eq. (fTSjl or ()14|) without considering the change in A with the field. A similar field 
dependence for a is also reported from another groupW . 

On the other hand, in KOS2O6 (Tc ~ 9.5 K), our preliminary /iSR data on a powder 
specimen exhibits a strong field dependence of A(/z)[nil, suggesting the presence of 
anisotropy in the order parameter. This is also consistent with the absence of coherence 
a peak in ^^K NMR-Ti^. 

5.5. Other examples 

As discussed in Section 3, quasiparticle excitation due to the Doppler shift is predicted 
to be stronger for a larger degree of manifoldness in the nodal structure of the order 
parameter. Thus, it is naturally expected that superconductors with (i-wave paring 
would exhibit a strong field dependence of \{h). This was proven to be the case by 
systematic /xSR studies on the vortex state of high-Tc cupratesjO]. A typical example 
is found in YBa2Cu306.95 (see Fig. El), in which 77 is reported to be 5-6.6 over a field 
range < if < 2 TjnSj. They later extended the measurement up to 7 T, where they 
found that the field dependence of A(/i) became weaker at higher fields (// ~ 2 for H > 2 
T) [21]. This is now understood to be a consequence of the non-local correction discussed 
earlier in Section 4. It must be noted that the non-local correction has a strong influence 
on the temperature dependence of A at higher fields. These results imply that one must 
be careful about the field range of measurements to evaluate the meaning of rj as well 
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Figure 9. Field dependence of the magnetic penetration depth (A) in YBa2Cu306.95 
(YBCO) and NbSe2 at T = O.SSTc, obtained by fitting with the modified London 
model, where Hc2 = 95 T and 2.9 T, respectively. A linear fit yields t] = 5-6.6 in 
YBCO and f .85(7) in NbSe2. A slightly smaller value {rj — 1.62(6)) is reported for the 
latter with an analysis by the GL model (after Ref.|3]). 



as the slope against temperature, d\/dT. 

The last example is NbSeo lMl lUj. where the situation is similar to that in YNi2B2C 
or MgB2. It has been suggested from the non-linear field dependence of 'y{h) that there 
must be excess quasiparticles induced by a magnetic fieldjSD]- However, the degree 
of non-linearity is considerably weaker than that observed in YNi2B2C, suggesting 
the smaller anisotropy in the order parameter. This has been supported by other 
experiments showing that NbSe2 has an anisotropic s-wave gap with Amax/Amin — 
2 jnSl inni EZI • As shown in Fig. 121 the absence of nodes, however, does not necessarily 
mean 77 = when the order parameter is anisotropic (or multi-gapped). It happens that 
a temperature of ~ 0.33Tc (where measurements were performed) IM] is relatively 
high, so that it is almost comparable to Amin- Thus, as discussed earlier, the region 
around Amin in the Fermi surface works effectively as nodes at such a high temperature. 
Given that this "quasi- node" scenario is correct, one can predict a larger non-linearity in 
7(/i) as well as a larger rj at higher temperatures. Such a tendency is actually reported 
in Ref . |94j . where they observed larger 77 at T = 0.6Tc. 

6. Summary and conclusion 

We demonstrated that the field dependence of the magnetic penetration depth \{h) 
provides a sensitive probe for quasiparticles induced by the Doppler shift. As 
summarized in Table 1, the slope 77 is positive when the superconducting order parameter 
has nodes (or a small gap equivalent to the node at a given temperature), while it is close 
to zero for the conventional isotropic order parameter. Despite the ambiguity associated 
with the slight dependence of A on the employed model for data analysis, the magnitude 
of rj provides a good measure for the degree of anisotropy. This would be useful in 
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selecting the paring symmetry and the associated mechanism of superconductivity for 
newly discovered materials. 

Table 1. Dimensionless parameter, rj, corresponding to the slope of X{h) against an 
external field obtained by fxSR. The column 'model' shows that of the field distribution 
used for each analysis: 'm-L' for the modified London model and 'G' for the Gaussian 
field distribution. T denotes the temperature where the field dependence of A was 
measured. The values for NbaSn are quoted from our preliminary report 98 . 



Tc (K) pairing symmetry rj model T (K) Hc2{T) (T) 



YNisBsC 


15.4 


anisotropic s (s + gl) 


0.95(1) 


m-L 


3.0 


7.0 


MgBs 


39 


double gap 


1.3(3) 


G 


10.0 


12.5 


NbSes 


7.0 


anisotropic s 


1.85(7) 


m-L 


2.3 


2.9 


YBa2Cu306.95 


93.2 


d 


5-6.6 


m-L 


31.0 


95 


CeRu2 


6.0 


isotropic s {h < 0.5) 


~ 


m-L 


2.0 


5.0 


Y(Nio.8Pto.2)2B2C 


12.1 


isotropic s 


~ 


m-L 


2.5 


4.0 


Cd2Re207 


1-2 


isotropic s 


^ 


G 


0.2 


0.37 


NbaSn 


18.3 


isotropic s 


~ 


G 


2.1 


~24 
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